Meson-meson scattering in lattice QED2+1 * 

H. R. Fiebig a and R. M. Woloshyn b 

a Physics Department, FIU University Park, Miami, Florida 33199, U.S.A. 

b TRIUMF, 4004 Wesbrook Mall, Vancouver, B.C., Canada V6T 2A3 

Scattering phase shifts of a meson-meson system in staggered 3-dimensional lattice QED are computed. The 
main task of the simulation is to obtain a discrete set of two-body energy levels. These are extracted from a 
4-point time correlation matrix and then used to obtain scattering phase shifts. The results for the 1 = and 
£ = 2 partial waves are consistent with short-range repulsion and intermediate-range attraction of the residual 
meson-meson interaction. 



1. INTRODUCTION 

The residual interaction between two hadron- 
like particles is of fundamental interest in strong- 
interaction physics. In a finite- volume lattice sim- 
ulation the interaction manifests itself in a dis- 
crete set of energy levels of the relative motion of 
the two particles. According to a recent proposal 
by Liischer [1] those energy levels can be used to 
extract scattering phase shifts. 

This program was implemented for a meson- 
meson system in lattice QED3. While this model 
is computationally much less demanding than lat- 
tice QCD4, it does have a confining phase and also 
exhibits chiral symmetry breaking. The shared 
phenomenology thus makes lattice QED3 a suit- 
able model for the problem at hand. 



2. LATTICE SIMULATION 

The sole purpose of the lattice simulation is to 
obtain the above-mentioned set of energy levels 
for the relative degrees of freedom. 

2.1. Lattice model 

We have used the compact version of the Wil- 
son action 
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where 0n denotes the plaquette angle. In three 
dimensions the staggered fermion action 
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describes two fermion species [2]. An external 
flavour index / = u,d was added to the Grass- 
mann fields. The usual staggered fermion matrix 
is called G~ l here. 

In order to match the phenomenology of QCD4 
it is essential that the field is in the confined 
phase. On a 24 2 x 32 lattice, which was used 
throughout the simulation, this is guaranteed at 
j3 = 1.5 with < W >= 0.03 being the value for 
the Wilson line. 

A second requirement is that the lattice be 
large enough to accommodate two mesons. To 
estimate the size of one meson, its mass m was 
obtained from simulations on L 2 x 32 lattices with 
various sizes L. The mass m(L) saturates expo- 
nentially at large L and is already within 5% of 
its asymptotic value at L = 10. 

A total of 64 gauge field configurations was gen- 
erated in the quenched approximation. 



2.2. Meson field operators 

A simple local one-meson operator is 

<j>p(t) = L- 2 J2^X d (x,t)Xu(x,t) 



(3) 



where p belongs to the set of L 2 discrete lattice 
momenta. Two mesons at center-of-mass momen- 
tum P = are then described by the operator 
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2.3. Time-correlation functions 

The threshold energy 2m of the noninteracting 
two-meson system is related to the square of the 
2-point correlation function 
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The interacting two-meson system is described by 
the 4-point correlation matrix 



C$(t,t ) 



(6) 



Both correlators C 1 - 2 - 1 and C 1 - 4 - 1 can be worked out 
in terms of contractions between the Grassmann 
fields contained in <f> and $. The simplistic flavour 
structure in (3) greatly reduces the number of 
terms in (6). There remain three basic diagrams 
that contribute to the 4-point correlation matrix 



C (4) = C (4A) + c 4B)_ 2ReC (4C) 

For example 
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This expression (and the two similar ones for 
(7( 4 ^) an d CJ^ ^) poses serious computational 
problems for two reasons. First, the entire prop- 
agator matrix G is needed. Second, the site sums 
have L 8 = 10 11 terms on the L = 24 lattice. They 
are numerically intractable. 

2.4. Random- source technique 

The matrix elements of C 1 - 4 - 1 can be computed 
efficiently with Fourier-modified random sources. 
Let Rg be complex Gaussian random vectors, 
then define 
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For each p and R the vector H can be obtained 
from a single call to, say, a conjugate-gradient 
routine. The 2-point and the 4-point correlators 



both have representations in terms of the auxil- 
iary vectors (9). For example 
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All contributions to C 1 - 4 " 1 require two independent 
random sources. The sums over i?' • > and R • ^ in 
(10) denote averages. For each average, 8 samples 
were used in the present simulation. 

2.5. Properties of the correlation matrix 

The eigenvalues X n (t,to) of the 4-point correla- 
tion matrix determine the two-body energy levels 
W n . We have [3] 



A n (Mo) = c n cosh[W n (t - t c )] 



(11) 



On the 24 2 x 32 lattice with periodic boundary 
conditions t c = 17. Figure 1 shows the time de- 
pendence of the five eigenvalues with the largest 
coefficients c n together with cosh fits according 
to (11). The matrix C 1 - 4 - 1 was truncated to mo- 
mentum indices p = (rii, n2)2ir/L with | «i 2 1 < 2. 
Certain symmetries of C 1 - 4 - 1 allow to restrict the 
range of momentum indices further to 13. 
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Figure 1. Eigenvalues of the 4-point correlator. 



Table 1 

Two-body strength factors and energy levels (Wo = 2m). 
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P 





0.958(25) 


1.032(08) 
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1.805(51) 


1.029(09) 
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0.191(06) 


1.114(09) 


0.047(3) 
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0.171(05) 


1.098(10) 


0.038(3) 
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0.043(01) 


1.209(12) 


0.117(6) 
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0.039(01) 


1.181(13) 


0.095(6) 



2.6. Two-body energy levels 

Table 1 shows the results of two-parameter fits 
to the eigenvalues according to (11). The thresh- 
old energy Wo = 2m was obtained from a fit to 
the square C^ 2 \t,to) 2 of the 2-point correlator. 

In principle, the relative energies 
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could now be compiled from the W n . In order 
to have better control over the errors, however, it 
is preferable to extract the k 2 ,n = 1...5, directly 
from fits to the ratios 



R„(t,to) = K(t,to)/C (2) (t,to) 2 



(13) 



The time dependence of the R n (t,to) is easily 
worked out using (11). With W n and c n fixed 
from the previous 2-parameter fits, the energies 
k 2 were obtained from 1-parameter fits to the ra- 
tios R n (t,to). The results are also shown in Ta- 
ble 1. The present data did not allow to make a 
fit to the n = 1 level. 

Comparison of Wo with W\ in Table 1 suggest 
the possibility that the n = 1 level may be a two- 
meson bound state. Although the present data 
are inconclusive in this regard, the residual in- 
teraction clearly seems to contain an attractive 
component. 

3. FINITE- VOLUME SCATTERING 

The discrete set of energy levels k 2 of the rel- 
ative motion, as listed in Table 1, is the basis 
for computing scattering phase shifts. Following 
Liischer [1] a nonlinear relation between the con- 
tinum phase shifts in partial wave £ = 0, 1...A 
can be derived by considering periodic replica- 
tions of the interaction potential and then match 



the physical (regular) solutions to the free (singu- 
lar periodic) solutions of the Helmholtz equation. 
Although the general outline of [1] remains un- 
changed, nearly all details have to be recalculated 
for the d = 2 dimensional case [4] . 

3.1. Cubic symmetry 

In d = 2 space dimensions the cubic symmetry 
group 0(2, Z) has 5 irreducible representations 
T = A\, A2, B\, B-2, E which are one dimensional 
with the exception of E, which has dimension 
two. The spaces £ spanned by Yi m (<f> 
with m = for £ = and m = ±1 for 
have the following decompositions 
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with E and E* being equivalent. 

3.2. Determinant condition 

The above-mentioned nonlinear condition for 
the phase shifts reads 



JJdet 



'2i6 



(r) 



M( r ) 



M(H + i 







(15) 



Here Af( r ) = M (v \q), with q = kL/2ir, is a ma- 
trix of size 2A + 1 (cut off at £ = A) that contains 
the scattering information of the cubic d = 2 geo- 
metry. The (diagonal) matrix e 2li comprises 
the phase shifts 6 £ ' (q) which are contained in 
the representation T, see (14). 

The task now consists in solving (15) for each 
of the relative energies q 2 , n = 2. ..5, listed in 



Table 1. Clearly the solution is not unique. We 
have treated each of the five factors in (15) sep- 
arately. Even then the solutions are not unique, 
in general. At an angular momentum cutoff of 
A = 8 the determinants have the sizes 3,2,2,2,4 
for T = A\, A2, B\, B2, E respectively. The so- 
lutions, however, become unique if one assumes 
that within each space T only one partial wave £ 

is dominant. This assumption [1] may be tested 

(r) 
by setting all but one 6\ to zero. In this way 
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one obtains 13 sets of phase shifts <5 
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4. THE PHASE SHIFTS 

We find that most of these phase shifts are in- 

(V) 

deed small, as expected [1]. In particular all 5\ 
with £ odd are essentially zero. Significant struc- 
ture is only found in the Ai-sector for £ = 
(s-wave) and the 52-sector for £ = 2 (d-wave). 
Those data are shown in Figures 2 and 3. 
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Figure 3. d-wave phase shifts (£=2) in sectors B^ 
(•) and 5i (o). 



In the d-wave the data have been plotted such 
that the phase decreases by 180° as k 2 becomes 
large. This is consistent with either a bound state 
(Levinson's theorem) or a sharp resonance close 
to the threshold. In either case, since the d-wave 
overlap tests somewhat larger distances (possibly 
the edges of the mesons), this feature would indi- 
cate intermediate-range attraction of the residual 
interaction. 

5. CONCLUSION 

The results of the meson-meson scattering sim- 
ulation in lattice QED3 possibly suggest short- 
range repulsion and intermediate-range attrac- 
tion of the residual force. These features allude 
to basic properties of hadronic interaction. 



Figure 2. s-wave phase shifts (£=0) in sector A\. 



The monotonic decrease of the s-wave phase 
shift, up to about k 2 = 0.1, is consistent with a 
repulsive residual interaction. Since the meson- 
meson overlap at short distances is largest in the 
s-wave, this feature may simply indicate Pauli re- 
pulsion. The data point at k 2 = 0.12 appears 
to indicate some other type of (possibly lattice) 
effect, which is not yet understood. 
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